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A method and apparatus for image processing 

r 

The present invention relates to methods and apparatus 
for image processing and, particularly, though not 
exclusively, to the" processing of Nuclear Magnetic 
Resonance image data. 

Certain imaging technologies, such as Magnetic Resonance 
(MR) imaging, permit the subject of the imaging process 
to be fully imaged in three spatial dimensions. They 
permit the acquisition of image data which is able to 
represent a 3D "view'' of the subject from any aspect or 
perspective in 3D space. This is in contrast to image 
data acquired through a series of predetermined 
"projections" which only enable "views' 7 limited to one of 
the predetermined projections employed in acquiring the 
image data. ; ^ 

Dynamic imaging, such as dynamic MR imaging, may involve 
the acquisition of a time-series of image data sets in 
which each data set is able to represent an image of the 
subject fully in 3D at one of a series of successive 
instants (or brief periods) in time. Dynamic imaging is 
often used to record the internal changes in* properties 
of a stationary subject which are induced" by a controlled 
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change-inducing influence. An example of this is 
contrast-enhanced medical MR imaging in which a "contrast 
agent" is introduced into a stationary subject (e.g. a 
person) which is detectable as an increase in the image 
5 contrast/brightness of those internal parts of the 
subject in which the agent is located. 

An analysis of the change in image properties (e.g. in 
image contrast/brightness), over time, of chosen fixed 

10 points or static regions within the imaged volume of the 
subject enables an assessment of the properties of the 
chosen fixed points/regions within that imaged volume. 
However, this assessment is rendered void if the chosen 
fixed points/regions of the imaged subject are assumed to 

15 be fixed/static but were, in fact, not static during the 
time-sequence in question. The result would be that 
changes in the properties of chosen points/regions would 
not be due to changes in the properties o*f the same 
point/region of the subject, but would instead be due to 

20 the. chosen point/region not representing the same 

point /region of the image subject for the duration of the 
time-sequence. Thus, image changes would be due to 
changes in the position of the subject (i.e. subject 
motion) rather than to changes in the subject per se. 



25 



For example, the image analysis procedure in breast 
dynamic imaging ' typically involves localisation of a 
lesion within an image of the subject containing a 
contrast agent. Image analysis may involve analysis of 
the shape of "time-intensity" curves . representing the 
temporal development of pixel intensity of chosen 
points/regions of the subject image and may involve 
pharmacokinetic modelling. All the techniques assume that 
the imaged volumes within the time-series are properly 
spatially aligned. A dynamic breast scan is a time 
consuming examination. The overall examination time can 
be up to 30 minutes. During such a long time a patient 
can move because of discomfort or tension. The motion 
most likely damages the spatial alignment of the imaged 
subject volumes within the time-sequence, and may lead to 
unreliable results, such as appearance of false lesions 
within the subtracted post-contrast image. ' Also, no 
reliable analysis based on- the ^time-intensity uptake 
curves is possible if the curves are not representative 
of the same point /region of the imaged subject. The 
motion of breast is not always easy to detect. Due to 
the high degree of elasticity of a breast tissue, a local 
motion may occur resulting in local deformation of some 
region the breast and tiny deformation of the remaining 
part of the volume. It follows that it is important to be 
able to assess a quality of alignment of image data, in 
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the local region of interest, as between imaged volumes 
within a time-sequence. 



The present invention aims to overcome at least some of 
5 the above deficiencies of the prior art. 

The present invention proposes to provide a method and/or 
apparatus enabling one . to determine if a 
contrast/intensity change of a chosen point within an 

10 imaged subject volume (e.g. a lesion) is a motion 
artefact. An apparatus and method for fast visual 
assessment of the imaged subject volume's alignment 
quality (within the time-sequence) is provided, and 
interactive display tool enabling effective use of this 

15 method is provided. 

In a first of its aspects, the present invention may 
provide a method of processing a time-sequence of 
separate image data sets each set comprising a plurality 
20 of image data items which each represent the location of 
an image pixel of the image subject according to a 
spatially fixed reference frame within which the subject 
is located, the method including the steps of: 

(a) selecting a plurality of separate image pixel 
25 locations within the reference frame; 

(b) selecting from each of a plurality of the 
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separate image data sets those image data items which 
represent image pixels located at the plurality of 
separate image pixel locations ; 

(c) generating a time-domain image according to the 
5 selected image data items wherein those image pixels of 
the time-domain image sharing a common image pixel 
location within the reference frame are arranged in 
-time-sequence in a line containing only those pixels. 

10 Thus, the reference frame within which the subject is 
located is common to all image data items of all image 
data sets of the time-sequence. By lining-up image 
pixels sharing the same spatial coordinates within the 
reference frame, one lines-up the temporally successive 

15 pixel representations of the same point of space (imaged 
at successive times within a time-sequence) within the 
reference frame in which the subject is located. In 
doing so, one does not necessarily line-up image pixels 
of the same points of the imaged subject since the 

20 subject (or a part of it) may have moved its location 
within the reference frame during time-sequence. 
However, should such subject movement occur, a 
consequence is often a visibly discontinuous change, or 
an irregularity, in a property of image pixels (e.g. 

25 brightness, or colour) as between neighbouring pixels of 
the line of pixels in the 'time-domain image. The 
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location of such a discontinuity/irregularity within the 
time-domain image identifies two image volumes within the 
time-sequence between which the subject motion occurred. 
This method may be applied to image data sets with image 
5 data items representing pixel locations extending only in 
2D (e.g. image "slices"), as well as those extending in 
all three spatial dimensions (i.e. image volumes). 

Preferably, each image pixel is an image voxel such that 
10 each image data item represents the location of an image 
voxel of the image subject according to a spatially fixed 
three-dimensional reference frame, preferably, each image 
data set contains image data items representing image 
voxel locations spanning all three dimensions of the 
15 reference frame. 

By selecting the .image data items representing a 
plurality of separate image pixel locations within the 
reference frame, a time-domain image comprising a 
20 plurality of separate lines each containing only those 

pixels sharing a common image pixel location may be 
•,._.^,.pr.ovidad,- ,.vJPhe,-ptorality of such separate lines are- 

preferably arranged side-by-side such that those image 
pixels of neighbouring lines which share a common image 
25 data set are also aligned so as to form a line. This 

arrangement assists in visually assessing image data for 
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the presence of motion artefacts since such artefacts are 
most likely to affect a plurality of image pixel 
locations within the reference frame rather than just one 
such location.. Accordingly, local image features may 
5 well display motion artefacts as a collective 

discontinuity or irregularity shared by several or each 
of the side-by-side lines within the time-domain image. 

The time-domain image is preferably generated such that 
10 • those image pixels of the time-domain image which are 
represented by image data items selected from a common 
image data set are arranged in a line containing only 
those pixels. 

15 Often small regions within the imaged subject volume 

suffer motion artefacts while the remaining volume does 
not. In order to facilitate easy visual identification 
of such localised motion, the plurality of separate image 
pixel locations selected from within the reference frame 

20 at step (aj are neighbouring such locations. 

Consequently, provided that localised motion affects a 
small region of spatda-1 dimensions -mat^in^^^axeeeding- 
the spatial separation between such neighbouring image 
pixel locations, then motion artefacts will appear in the 

25 time-domain image (generated at step (c) ) as a collective 
discontinuity/irregularity affecting the same regions of 
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neighbouring of the side-by-side lines therein. Such 
collective discontinuity not only indicates the extent of 
subject motion but is also easier to spot within a time- 
domain image. 

5 

The time-domain image preferably comprises a matrix of 
pixels arranged to define columns and rows, wherein each 
column (or, alternatively, each row) contains only those 
image pixels sharing a common pixel location within the 
10 reference frame, and each row (or, alternatively, each 

column) contains only those image pixels selected from a 
common image data set. 




Preferably, the time-domain image is generated such that 
15 each separate line of image pixels of a common image 
pixel location contains successive line segments 
containing only those pixels which share a common image 
data set. Thus, each such line segment is comprised only 
of several copies of the same image pixel. Each segment 
20 within a line may contain the same number of image pixels 
(copies) thereby rendering each line segment the same 

.-Length*«*Motion artefacts in a single image .pixel.-,..-. 

location result in a visible shift of an entire line 
segment. This is much easier to identify than is a shift 
25 in a single pixel. More preferably, each of the side-by- 
side lines of the time-domain image comprises successive 
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such line segments, and preferably corresponding line 
segments of each such side-by-side line are equal in 
length (i.e. the same number of image pixels therein). 

5 Preferably, the selection of the plurality of image pixel 
locations at step (a) comprises selecting a target image 
pixel location within the reference frame and 
subsequently selecting all image pixel locations which 
are intersected by. a predefined first locus . intersecting . 

10 the target image pixel location. The locus may be any 
predefined shape and may be defined by a mathematical 
expression in terms of the coordinates of the reference 
frame, or may simply be a preselected trail of successive 
image pixel locations extending through the reference 

15 frame. 

Preferably, the predefined first locus is a line 
extending through the reference frame. The line 
preferably extends through only one dimension of the 
20 reference frame and intersects all pixel locations within 
a column or a row of pixel locations within the reference 
. frame. In. this way„ : a. IjLnear cross-sectipna-l view (in one 
dimension, i.e. ID) through the imaged subject volume is 
provided. 
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For example, method may include: selecting a target one 
of the plurality of separate image data sets and 
selecting therefrom those image data items which 
represent image pixels located in a common first plane 
5 containing the first locus; and generating a first 
spatial-domain image according to the data items so 
selected. The method preferably also includes 
generating a composite image containing the first 
spatial-domain image and the time-domain image aligned 
10 such that those pixels of the first spatial-domain image 
and of the time-domain image which share a common image 
pixel location within the reference frame are arranged 
in-line. • 

15 Such in-line arrangement enables a view of a spatial 
cross-section of the imaged subject volume and a 
simultaneous view of a time-domain image of the subject, 
each view sharing at least one common spatial dimension 
defined, by the first locus. The time-domain image also 

20 extends in the time dimension while the first spatial- 

domain image extends in two spatial dimensions. In-line 
.-arrangement of- the two . images • enables a direct- vi*»l*. 
comparison of the image features present along the common 
spatial dimension shared by the two images. 
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Preferably the first locus is simultaneously visually 
presented within the composite image as a first locus 
image overlaying one or each of the spatial-domain and 
the time-domain images and intersecting those image 
5 pixels of the respective images which represent pixel 
locations within the reference frame intersected by the 
first locus. Alternatively, such image pixels may be 
highlighted or otherwise rendered distinct from their 
neighbouring image pixels within the composite image. 

10 

Preferably, the selection of the plurality of image pixel 
locations at step (a) comprises selecting the target 
image pixel location within the reference frame and 
subsequently selecting all image pixel locations which 

15 are intersected by a predefined second locus intersecting 
, the first locus at the target image pixel location. The 
second locus may be any predefined shape and may be 
defined by a mathematical expression in terms of the 
coordinates of the reference frame, or may simply be a 

20 preselected trail of successive image pixel locations 
extending through the reference frame. 

Consequently, the composite image may contain a second 
time-domain image generated according to any one of the 
25 method steps discussed above in respect of the second 

predefined locus which intersects the first locus. The 
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second locus is preferably, substantially coplanar with 
the first locus. The first spatial-domain image and the 
second time-domain image are preferably aligned such that 
those pixels of the first spatial-domain image and of the 

5 second time-domain image which share a common image pixel 
location within the reference frame are arranged in-line. 
Preferably, a second locus image of the form of the first 
locus image may simultaneously be presented in either or 
each of the spatial-domain image and the second time- 

10 domain image. 

The first predefined locus is preferably a line and the 
second predefined locus is preferably also a line 
substantially perpendicular to the first predefined 
15 locus. This provides a cross-hair effect in a locus 

image comprising the first and second locus images. The 
cross-hairs will centre upon the target pixel location 
within the spatial-domain image. 

20 For example, when the first and second locus' are 
perpendicular and coplanar the composite image, 
comprising the first spatial-domain image -and-uhoth:., the- - 
first and second time-domain images, provides a composite 
image representing three dimensions in which each time- 

25 domain image shares a separate one of the two spatial 
dimensions of the first spatial-domain image, and in 
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which each time -domain image extends in the time 
dimension. Graphically, the time dimension of each time- 
domain image may extend perpendicularly to the spatial 
dimension of the image in question. A further spatial 
5 dimension orthogonal to those represented in the first - 
spatial-domain image may also be represented within the 
composite image thereby rendering the latter 
representative of four dimensions (i.e. three spatial and 
one temporal) . 

10 

Thus, preferably, the selection of the plurality of image 
pixel locations at step (a) comprises selecting the 
target image pixel location within the reference frame 
and subsequently selecting all image pixel locations 

15 which are intersected by a predefined third locus 

intersecting the first locus and the second locus at the 
target image pixel location. The third locus may be any 
predefined shape and may be defined by a mathematical 
expression in terms of the coordinates of the reference 

20 frame, or may simply be a preselected trail of successive 
image pixel locations extending through the reference 
v * 'f ranie i i , cwwt,.,i-* l • ^^at^.^^-:^^' " * ** ' * : " 

Consequently, the composite image may contain a third 
25 time-domain image generated according to any one of the 
method steps discussed above in respect of the third 
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predefined locus which intersects the first and second 
locus'. The third locus is preferably non-coplanar (e.g. 
perpendicular) with the first locus and the second locus. 

5 The method may further include: selecting from the target 
one of the plurality of separate image data sets those 
data items which represent image pixels located in a 
common second plane being non-parallel and non-coplanar 
with the first plane and containing the first locus; 

10 generating a second spatial domain image according to the 
data items so selected. The method preferably further 
includes generating the composite image so as to contain 
the second spatial-domain image aligned such that those 
pixels of the second spatial-domain image and of the 

15 time-domain image which share a common pixel location 
within the reference frame are arranged in-line. 

In this way, a third spatial domain is represented within 
the composite image by the presence of the second 

20 spatial-domain image which also shares one common spatial 
dimension with both the first spatiai-domain image and 
the .first time-domain •.iiaag.e.--._Ih±s., -common. .spa&La-l 
dimension is defined by the first locus and preferably 
the common spatial axis of each of the first time-domain, 

25 first spatial-domain and second spatial-domain images are 
aligned in parallel and side-by-side such that the time 
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axis of the first time-domain image and the spatial axes 
of the first and second spatial-domain images which are 
perpendicular to the first locus, are all aligned 
sequentially in-line and in parallel. This permits rapid 
5 " visual assessment of image data in three spatial and one 
temporal dimension. 

The method may also include: selecting from the target 
one of the plurality of separate image data sets those 

10 data items which represent image pixels located in a 

common third plane being non-parallel and non-coplanar 
with the first plane (and the second plane, preferably) 
and containing the second locus; and generating a third 
spatial domain image according to the data items so 

15 selected. Preferably the method 6 further includes 

generating the composite ittiage so as to contain the third 
spatial-domain image aligned such that those pixels of 
the third spatial-domain image and of the second time- 
domain image which share a common pixel location within 

20 the reference frame are arranged in-line. Again, most 

preferably the parallel alignment of common spatial axes 
- ^..^tis^praasM^ arrangement as- discus-sed — -v< - ^~./r* * 

above - 

25 The first, second and third planes are preferably all 
mutually substantially perpendicular. Preferably, the 



16 

second plane and the third plane, in respect of which the 
second and third spatial-domain images are generated, 
each contain the third locus. 

5 The third time-domain image and either the third or the 
second spatial-domain image are preferably aligned such 
that those pixels of the third time-domain image and of 
either the second or the third spatial-domain image which 
share a common image pixel location within the reference 

10 frame are arranged in-line'. The spatial axis of the 

third time-domain image may be arranged parallel and in- 
line with the spatial axis of the first time-domain image 
(i.e. the two images side-by-side) when the third time- 
domain image is aligned with the third spatial-domain 

15 image. Alternatively, the spatial axis of the third time- 
domain image may be arranged parallel and in-line with 
the spatial axis of the second time-domain image (i.e. 
the two images side-by-side) when the third time-domain' 
image is aligned with the second spatial-domain image. 

20 Preferably, a third locus image of the form of the first 
or second locus image may simultaneously be presented in 
either. oa?_~ each- ©f v Xh^aecond^ 
image and the third time-domain image. 

25 Preferably, all image pixels within each of the images 
contained within the composite image which are 
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intersected by a locus are highlighted so as to render 
them distinguishable from their neighbouring image 
pixels, or are overlaid by a locus image which represents 
the position of a locus within an image . 

Planes are non-parallel if the normal extending from one 
plane is not parallel to the normal extending from 
another "non-parallel" plane. 

In a second of its aspects the present invention may 
provide apparatus for processing a time-sequence of 
separate image data sets each set comprising a plurality 
of image data items which each represent the location of 
an image pixel of the image subject according to a 
spatially fixed reference frame within which the subject 
is located, the method including the steps of: 

(a) first selection means selecting a plurality of 
separate image pixel locations within the reference 
frame ; 

(b) second selection means for selecting from each 
of a plurality of the separate image data sets those 

:Aima.ge^-da-ta£i i.t-ems which represent image pixels* loeafegcU-afc 
the plurality of separate image pixel locations; 

(c) image generating means for generating a time- 
domain image according to the selected image data items 
wherein- those image pixels of the time-domain image 
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sharing a common image pixel location within the 
reference frame are arranged in the time-sequence in a 
line containing only those pixels. 

5 Each image pixel is preferably an image voxel such that 
each image data item represent the location of an image 
voxel of the image subject according to a spatially fixed 
three-dimensional reference frame. Preferably each image 
data set contains image data items representing image 

10 voxel locations spanning all three dimensions of the 
reference frame. 

The first selection means is preferably arranged to 
select a target image pixel location within the reference 
15 frame, and includes locus means for defining a first 
locus intersecting the target image pixel location, 
wherein the first selection means is arranged to select 
all image pixel locations which are intersected by the 
first locus. 

20 

The first locus is preferably a line extending through 
. --.the easencaaL JEraras«- •• Saa^ine -preferably extends through 
only one dimension of the reference frame and intersects 
all pixel locations within a column or a row of pixel 
25 locations within the reference frame. 
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The image generating means is. preferably arranged to 
generate the time-domain image such that those image 
pixels of the time-domain image which are represented by 
image data items selected from a common image data set 
are arranged in a line containing only those pixels. 

The image generating means is preferably arranged to 
.generate the time-domain image comprising a matrix of 
pixels arranged to define columns and rows, wherein each 
column (or, alternatively, each row) contains only those 
image pixels sharing a common pixel location within the 
reference frame, and each row (or, alternatively, each 
column) contains only those image pixels selected from a 
common image data set. 

Preferably -the second selection means is arranged to 
select from a target one of the plurality of separate 
image data sets those image data items which represent 
image pixels located in a common first plane containing 
the first locus; wherein the image generating means is 
arranged to generate a first spatial-domain image 
according to the data items so., selected*; Er.ef erab±^t.h^.rt 
image generating means is further arranged to generate a 
composite image containing the first spatial-domain image 
and the time-domain • image aligned such that those pixels 
of the first spatial-domain image and of the time-domain 
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image which share a common image pixel location within 
the reference frame are arranged in-line. 

The locus means is preferably arranged to define a second 
5 locus which intersects the first locus at the target 
image pixel location and is preferably substantially 
coplanar with the first locus. The second selection means 
is preferably arranged to select all image data items the 
pixel locations of which are intersected by the second 
10 locus. The image generation means is preferably arranged 
to generate a second time-domain image in respect of 
image data it.ems so selected. Preferably the image 
generating means is further arranged to generate a 
composite image which contains a second time-domain image 
15 in which the first spatial-domain image and the second 
time-domain image are aligned such that those pixels of 
the first spatial-domain image and of the second time- 
domain image which share a common image pixel location 
within the reference frame are arranged in-line. 

20 

The locus means is preferably arranged to define the 
firsts.-, locus --as^iar-iritne.. and the -second locus is a. line. - 
substantially perpendicular to the first locus. 



25 The locus means is preferably arranged to define a third 
locus which intersects the first locus and the second 
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locus at the target image pixel location and is contained 
within one of the second plane and the third planes. The 
second^ selection means is preferably arranged to select 
all image data items the pixel locations of which are 
intersected by the third locus. The image generation 
means is preferably arranged to generate a third time- 
domain image according to image data items so selected. 
Preferably the image generating means is further arranged 
. to generate a composite image 
which contains the third time-domain image, wherein the 
third time-domain image and one of the second and the 
third spatial-domain images respectively are aligned such 
that those pixels of the third time-domain image and 
those pixels of the one of the second and the third 
spatial-domain images which share a common image pixel 
location within the reference frame are arranged in-line. 

The second selection means is preferably arranged to 
select from the target one of the plurality of separate 
image data sets those data items which represent image 
pixels located in a common second plane being non- 
parallel: and non-coplanar- with t;he._f irsfcf;.plane -.aaekrv^ • . 
containing the first locus; the image generating means is. 
preferably arranged to generate a second spatial-domain 
.image according to the data items so selected. 
Preferably the image * generating means is also arranged to 
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generate the composite image so as to contain the second 
spatial-domain image aligned such that those pixels of 
the second spatial-domain image and of the time-domain 
image which share a common pixel location within the 
5 reference frame are arranged in-line. 

The second selection means, is preferably arranged to 
select from the target one of the plurality of separate 
image . data sets those data items which represent image 

10 pixels located in a common third plane being non-parallel 
and non-coplanar with the first plane (and preferably the 
second plane) and containing the second locus; and the 
- image generating means is preferably arranged to generate 
a third spatial domain image according to the data items 

15 so selected. Preferably the image generating means is 
also arranged to generate the composite image so as to 
contain the third spatial-domain image aligned such that 
those pixels of the third spatial-domain image and of the 
second time-domain image which share a common pixel 

20 location within the reference frame are arranged in-line. 

«>» .^-:>^The^±rnfager.:g-enerating means . is pref erably- operable- -to- 

generate a composite image, in which the spatial axis of 
the third time-domain image may be arranged parallel and 
25 in-line with the spatial axis of the first time-domain 
image (i.e. the two images side-by-side) when the third 
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time-domain image is aligned with the third spatial- 
domain image. Alternatively, the spatial axis of the 
third time-domain image may be arranged parallel and in- 
line with the spatial axis of the second time-domain 
5 image (i.e. the two images side-by-side) when the third 
time-domain image is aligned with the second spatial- 
domain image. Preferably, the locus means is operable to 
define a third locus image of the form of the first or 
second locus image and to simultaneously present the 
10 third locus image in either or each of the second or 
third spatial-domain image and the third time-domain 
image . 

The second selection means is preferably arranged to 
15 select the first, second and third planes to be mutually 
substantially perpendicular. 

A plurality of items are "in-line" if they are all 
arranged in/along a common single straight line. 

20 

The invention may also comprise computer means programmed 
to perform the method according :to- the invention is its 
second aspect including some, none or all of the 
preferred features or variants discussed above. 
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In a third of its aspects, the present invention may 
provide computer means programmed to perform the method 
according to the invention in its first aspect, including 
none or all of the preferred features or variants 
5 discussed above. 

In a fourth of its aspects, the present invention may 
provide a computer program product containing a 
computer program for performing the method according to 
10 the invention in its first aspect, including none or all 
of the preferred features or variants discussed above. 

The present invention may also provide a computer program 
containing computer code which, when implemented, 
15 performs the method according to the invention in its 
first aspect, including none or all of the preferred 
features or variants discussed above. 

The present invention may also provide an image generated 
20 according to the invention in its first aspect, including 
none or all of the preferred features or variants 
•discussed above. : : - .*•".* 

The invention shall now be described in terms of the 
25 following non-limiting examples with reference to the 
accompanying drawings in which: 
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Fig. 1 illustrates a plurality (a) to (i) of time- 
domain images displaying various degrees of image 
discontinuity/irregularity due to subject motion (i.e. 
motion artefacts) ; 
5 Fig. 2 illustrates two time-domain images (a) and 

(b) displaying the temporal development of contrast agent 
uptake within an imaged subject at selected regions 
thereof; 

Fig. 3 illustrates a composite image -comprising 
10 three time-domain images and three spatial-domain images 
defined by three perpendicular locus'; 

Fig. 4 illustrates schematically a reference frame 
within which a subject is located and according to which 
the image data items of the imaged subject represent the 
15 locations of image voxels representing points within the 
space defined by the reference frame; 

Fig. 5 illustrates schematically a motion artefact 
within a time-domain image generated according to a 
predefined locus L. 

20 

Referring to Fig. 4 there is illustrated a reference 
.frame R comprising t.h.ree... orthogonal axes (x, y, z) 
defining a spatially fixed coordinate system within which 
is located a subject O in respect of which image data 
25 items fully representing the object in 3D have been' 
acquired. The acquired 3D image data items each 
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represent the location of an image voxel of the image 

subject 0 according to the spatially fixed reference 

frame R within which the subject is located. A data set 

» 

comprising a plurality of such image data items 
5 accordingly represents a plurality of locations within 
the reference frame R spanning a predefined volume of 
space V. A given such data set is apquired over a 
predetermined brief period of time sufficient to acquire 
.all of the image data items comprising the set. 
10 Acquisition of such image data sets may be repeated a 
predetermined number of times (e.g. M times) thereby 
defining a time-sequence of separate image data sets each 
representing the same fixed volume V at the same fixed 
image voxel locations within the reference frame R. 

15 

Given such a time-sequence or series of M dynamic data 
sets (volumes) , each imaged volume defined by a given 
data set within the sequence is represented as a stack of 
frames with each frame comprised of a 2D lattice of image 
20 voxels. The resulting image volume comprises a lattice 
of image voxels having dimensions N x xN y xN z (respectively 
number of columns and rows in a frame,- "and a stack 
length) . 
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The method includes the initial step of selecting a 
plurality of separate image voxel locations within the 
reference frame (R) by: selecting a target voxel from the 
volume lattice V (hereinafter referred to as the "seed 
5 voxel") ; selecting from the volume V a frame (P) 

containing the seed voxel; and, from the selected frame 
selecting the row (L) containing the seed voxel. The 
selected row represents a predetermined linear locus with 
which to select image voxel locations (once the seed 
10 voxel is identified) . The locus need not be linear. 

By selecting from each of the M separate image data sets 
those image data items which represent image voxels 
located at the plurality pf selected image voxel 

15 locations defined by the locus L, one obtains image data 
items for use in generating a time-domain image. This 
data selection step may occur subsequent to selection of 
the image voxels locations via locus L, or may occur in 
tandem by, for example, after selecting a seed voxel, 

20 subsequently selecting from each of the M data sets the 

image data items contained in the selected frame (P) , and 
.then selecting. image .data items intersected by the. lpcus . 
L within each selected frame of data. 

25 The process of selecting voxel locations may be 

undertaken with a computer and visual display unit 
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implementing software according to the invention enabling 
a seed voxel (S) , a frame (P) and a row or other locus 
(L) to be selected interactively from a visual display of 
an image or images of the object derived from one, some 
5 or all of the M image data sets. 



Subsequent to data selection, a time-domain image is 
generated from the selected image data items wherein 
those image voxels of the time-domain image sharing a 
10 common image voxel location within the reference frame 
are arranged according to the time-sequence in a line 
containing only those voxels. 

This may be done by combining the selected rows (R) into 
15 an image where order of rows follows the order of 
acquisition of the M data sets (volumes) from which the 
rows of data were selected. The resulting time-domain 
image has dimensions N x xM . 

20 Any point on the time-domain image corresponds to a point 
in the reference frame R as represented by a image data 
set representing the- volume V- - Given, the coordinates of 
the seed voxel S are (x s> y S9 z s ) , the points P TD (u,v) within 
the volume V which are represented by the time-domain 
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image corresponds to the four-dimensional points in the 
original data: 

Pad (*> y,*>t) = Pad ( w > J>5 > z s > v ) ( 1 > 

5 

Where x,^^ are volume lattice coordinates in the 
reference frame R, and the fourth coordinate, t, is a 
temporal coordinate, i.e. index of the dynamic volume 
within the time-sequence. 

10 

Using the same principle, one may create a time-domain 
image selecting columns L2 passing through the seed voxel 
S instead of the rows L. That is. to say, the locus 
according to which image voxel locations are selected may 
15 be a column L2 rather than (or in addition to) a row L. 
In such a case the points of interest represented by the 
time-domain image will be defined by: 

Pad ( x > y> *) = Pad( x s> v > z s> u ) ( 2 > 

20 

Furthermore, one may create a time-domain image selecting 
a ^de^th 7 '** 5 column L3 passing through the seed voxel S' 
perpendicular to the row L and column L2 . That is to say, 
the locus according to which image voxel locations are 
25 selected may be a column L3 rather (or in addition to) 
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than a row L, or column L2 . In such a case the points of 
interest represented by the time-domain image will be 
defined by: 

5 P AD (x>y,z 9 t) = P 4D (x S9 y s9 u ,v) (3) 

This equation assumes that pixels selected from a common 
data set are in rows. The choice of image plane from 
which this third locus L3 is selected may be the plane 
10 parallel to the plane containing the z axis and the x 
axis of the reference frame R, or may be the plane 
perpendicular to that and containing the z axis and the y 
axis of the reference frame. 



15 Figure 5 schematically illustrates a motion artefact 
within a time-domain image generated according to a 
predefined locus L. The time-domain image comprises a 
plurality of separate vertical lines 1,2 and 3 etc. 
(image columns) each containing only those voxels sharing 

20 a common image voxel location. The plurality of such 
separate vertical lines are arranged side-by-side such 
that those image . voxels -sharing a common image data set 
aligned so as to form a horizontal line (image rows) . 
This arrangement assists in visually assessing image data 

25 for the presence of motion artefacts since such artefacts 
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are most likely to affect a plurality of image voxel 
locations within the reference frame rather than just one 
such location. Accordingly, local image features may , 
well display motion artefacts as a collective 
discontinuity or irregularity shared by several or each 
of the side-by-side vertical lines 1, 2, 3 etc. (image 
columns) within the time-domain image. 

The time-domain image of Figure 5 is generated such that 
those image voxels of the time-domain image which are 
represented by image data items selected from a common 
image data set (e.g.- set number m or set number m+1 
within the time-sequence of M sets) are arranged in image 
rows containing only those voxels. Consequently, the 
time-domain image comprises a matrix of voxels arranged 
to define columns and rows, wherein each column contains 
only those image voxels sharing a common voxel location 
within the reference frame R, and each row contains only 
those image voxels selected from a common image data set 
(e.g. set number m) . 

In "order to facilitate easy visual identification of 
localised subject motion, the plurality of separate image 
voxel locations selected from within the reference frame 
R are neighbouring such locations. Consequently, 
provided that localised motion affects a small region of 
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spatial dimensions matching or exceeding the spatial . 
separation between such neighbouring image voxel 
locations, then motion artefacts will appear in the time- 
domain image as a collective discontinuity/irregularity 
5 affecting the same regions of neighbouring image columns 
l f 2 and 3 etc. Such collective discontinuity not only 
indicates the extent of subject motion but is also easier 
to spot within a time-domain image. 

Each separate image column 1,2,3 etc. of voxels contains 
10 successive column segments (e.g. segment 1A and segment 
IB of column 1) containing only those voxels which share 
a common image data set. Each such line segment is 
comprised only of 3 copies of the same image voxel (any 
number of copies could be used in practice) . Motion 

15 artefacts in a single image voxel location result in a 
visible shift of an entire line segment. This is much 
easier to identify than is a shift in a single voxel.' 
Thus, while segment 1A of image column 1, derived from 
the data set number m displays an image feature 

20 represented by triangles in Figure 5, segment lb of image 
column 1, derived from the data set number m+1 displays 
an image f eatune^-xepresented by circles. However, the 
segment of column 2 associated with set m also displays 
an image feature represented by circles. This indicates 

25 that, in between acquiring image data set m and set m+1, 
the object feature responsible for the image feature 
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represented by circles has moved by an amount equal to 
the separation between the image voxel locations 
represented by columns 1 and 2 of the time-domain image. 
Thus, the visible discontinuity in columns 1 and 2 (and 
3) of the time-domain image is a motion artefact and not 
a property of the object per se. 

Figures 1(a) to l(i) illustrate a series of time-domain 
images of a breast imaged according to the present 
invention with data acquired using MR. Figure 1 shows 9 
separate examples (Figs. 1(a) to l(i)) of horizontal time- 
domain images each for a time-sequence comprising seven 
separate breast measurements {i.e. M=7) . The dynamic 
series is derived from seven image volumes: two pre- and 
five post-contrast , i.e. image data acquisition performed 
before (twice) and after (5 times) administration of the 
contrast agent Gd-DTPA. The time interval between 

acquisition of successive data sets was 90 sec. The data 
sets were acquired using 3D Tl-weighted fast spoiled 
gradient echo sequence. For better visual representation 
the columns of the time-domain images generated with 
column segments each containing two. identical, image 
voxels (see figure 5 and the discussion thereof) . 

Looking at the examples, one can easily conclude that 
cases (a-c) show no motion, in cases (d-g) motion is 
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obvious, and cases (h, i) present local motion (circled 
on the images) . 

The time-domain imaging approach may be used for 4- 
5 Dimensional visualisation of contrast enhancement of a 
subject containing a contrast agent. This application of 
the time-domain image preferably employs "subtracted" 
images in which the voxel value of each "subtracted" 
image data item represents the difference in voxel 
10 intensity (i.e. the contrast) as between the same voxel 
position in two different images, the first image being 
of the subject containing no contract agent and the 
second image being of the subject after administration of 
a contrast agent. In this case a time-domain image 
15 created for a seed voxel S selected inside a lesion on 
one a subtracted image will show the temporal pattern of 
contrast enhancement . 

Figure 2 shows fragments of horizontal time-domain MR 
20 breast images. The images display two lesions (Fig. 2 (a) 
and Fig. 2(b)) with different patterns of voxel 
enhancement. In both case the seed voxel (S) was selected 
near the centre of a lesion. The time axis is oriented 
from top to the bottom of the images. This means that the 
25 top row of each image was selected from the first post- 
contrast image data set, • the next one- from- "the second, 
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and so on. One can see in the figure 2(a) that the 
central part of the lesion is not enhanced on the first 
image (top line of figure 2(a) - inticated by arrow), but 
becomes enhanced at the bottom of the image (i.e. at a 
5 later time within the time-sequence) . At the same time, 
the peripheral region on the lesion is enhanced at the 
• top of the image of figure 2 (a) . The figure- 2 (b) presents 
opposite behaviour when maximum of the lesion central 
region enhancement is at the top, and maximum of the 
10 periphery enhancement is at the bottom, when the contrast 
is already washed out from the centre (pointed to by an 
arrow in fig, 2(b)). 

Time-domain images may be generated which represent the 
15 time-development of image voxels intersected by 
perpendicular linear locus 9 lying within the same 
selected frame/plane (P) . For example the time-domain 
image derived from the row locus (L) may represent the 
horizontal dimension of the volume V, while the image 
20 derived from the column locus (L2) represents the 
vertical dimension. A third complimentary time-domain 
image may be created - using similar approach according to 
a third locus L3 extending through the seed voxel S in a 
direction orthogonal to the frame/plane (P) , and would 
25 represent the depth dimension. 
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Equations 1, 2 and 3 imply a way of 4D navigation of a 
dynamic series of images. A 4-Dimensional orthogonal 
display and an interactive display tool for 4-D 
visualisation, navigation and alignment quality control 

5 of the time-sequence of 3D image data sets is provided 
according to the present invention. The tool uses the 
orthogonal display to select a point of interest - the 
seed voxel/voxel- using one of the M image data sets 
volumes at a time. Once a seed voxel S is .chosen, 

10 horizontal, vertical and depth time-domain images are 
created and displayed as discussed above. The selection 
of a new point of interest (i.e. a new seed S) on any one 
of the time-domain images results in the reselection, 
from the time-series of image data sets, of the 

15 appropriate image data set to be used for a concurrent 
display of the image plane (P) containing the reselected 
seed voxel, along with a consistent change in the spatial 
point of focus (i.e. the seed) within the image. 

20 Figure 3 illustrates a composite interactive image 

display 30 comprising a first spatial-domain image 31 
generated from image data items . selected from a .single 
image data set (e.g. set number m) which represent image 
voxels located in a common first plane (P) containing a 

25 first locus 45. The interactive display generates a 

composite image containing the first spatial-domain image 
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31 and a horizontal time-domain image 33 aligned such 
that those voxels of the first spatial-domain image and 
of the time-domain 33 image which share a common image 
voxel location within the reference frame are arranged 
5 in-line. 

The first locus is simultaneously visually presented 
within the composite image as a first locus image 
overlaying each of the spatial-domain 31 (item 45) and 
10 the time-domain 33 images (item 4 3) and intersecting 
those image voxels* of the respective images which 
intersect the first locus in the reference frame. 

The composite image contains a vertical time-domain image 
15 35 in respect of a second locus 44 which intersects the 
first locus 45 at the seed voxel, and is substantially 
coplanar with the first locus. The first spatial-domain 
image 31 and the second time-domain image 35 are aligned 
such that those voxels of the first spatial-domain image 
20 and of the second time-domain image which share a common 
image voxel location within the reference frame are 
arranged in-line. A second locus image (44 , 39) of the 
form of the first locus image is simultaneously presented 
in each of the spatial-domain image 31 and the second 
25 time-domain image 35. 
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The first locus is an image row within the first spatial- 
domain image 31 and the second locus is an image column 
perpendicular to the first predefined locus. This 
provides a cross-hair effect in a locus image comprising 
5 the first and second locus images. The cross-hairs 
centre upon the seed voxel location within the first 
spatial-domain image 31. 

The composite image 30 further includes a second 32 and a 
10 third 36 spatial-domain image each separately generated 
from image data items selected from the same image data 
set employed for the first spatial-domain image. The 
image data items employed in generating the second 
spatial-domain image represent image voxels located in a 
15 common second plane which contains the first locus 4 5 
(see item 37 of image 32) . The image data items 
employed in generating the third spatial-domain image 
represent image voxels located in a common third plane 
which contains the second locus (see items 38 and 36) . 
20 Each of the second and third planes is non-parallel and 
non-coplanar with the first plane P and with each other. 

The composite image contains the second spatial-domain 
image 32 and the third spatial-domain image 36 each 
25 aligned such that those voxels of the first spatial- 
domain image, the second spatial-domain image and the 
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first time-domain image which share a common voxel 
location in the reference frame are arranged in-line. 
Similarly, those voxels of the first spatial-domain 
image, the third spatial-domain image and the second 
5 time-domain image which share a common voxel location 
within the reference frame are arranged in-line. 
Consequently, a third spatial domain (i.e. depth) is 
represented within the composite image by the presence of 
the second and/or third spatial-domain images. The 

10 second spatial-domain image shares one common spatial 

dimension (i.e. the horizontal dimension) with both the 
first spatial-domain image 31 and the first time-domain 
image 33. Similarly, the third spatial-domain image 
shares one common spatial dimension (i.e. the vertical 

15 dimension) with both the first spatial-domain image 31 

and the second time-domain image 35. A third time-domain 
image 34 is provided within the composite image 30 
aligned to the horizontal axis of the third spatial- 
domain image 36. 

20 

Each image within the composite image contains a locus 
image (items. 37,. 38, 39, 40, 41, 42, 43, 44, 45) 
representing those voxels within the image which are 
intersected by a locus. 



25 
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Figure 3 presents the general view of a 4-Dimensional 
orthogonal display according to the present invention. 
The display shows three planes in orthogonal reformat of 
a coronal stack of subtracted post-contrast breast 
5 images. A vertical time-domain image 35 is displayed next 
to the sagittal frame 36, and a horizontal time-domain 
image 33 - below the axial frame 32. The "depth" time- 
domain image 34 is displayed at the bottom right corner 
of the composite image 30, and has a horizontal axis, and 
10 locus image 40, aligned with that of the sagittal frame 
36 and locus image 38. Similarly, the vertical axis of 
the "depth" time-domaim image 34 is aligned with that of 
the horizontal time domain image 33. 

15 The relative positions of the separate images within the 
composite image 30 are selected in such a way that 
aligned axes coincide in the 4D space. The image on 
display in Figure 3 is derived from an image data set 
representing the first post-contrast . image of the 

20 subject. The display is focused on the central region of 
the lesion. The time-domain images show a focal pattern 
of enhancement (i.e. central part of the lesion that is 
not enhanced on the first post-contrast image, but 
becomes enhanced later.) 
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A subtraction image is usually employed in dynamic 
imaging of breast tissues for localisation of a lesion, 
and visual analysis of lesion morphology. It follows that 
5 the quality of the subtraction image is essential for the 
correct interpretation of the image data. 

The present invention is effective for patient motion 
detection which, if left undetected, could lead to poor 
10 or erroneous image diagnosis. 

The present invention provides a method and apparatus for 
the assessment of the alignment of 3D imaged subject 
volumes and is suitable for detecting motion artefacts 

15 resulting from local areas of motion/ deformations within 
the imaged volume as well as global motion. It is of 
particular importance for breast imaging because of the 
high elasticity of breast tissue which renders it 
susceptible to localised (non-global) movements due to 

20 deformations . 

....... . Furthermore, the present invention also provides a method 

of visualisation of dynamic contrast enhancement in a 
subject. A time-domain image of subtracted post-contrast 
25 images (or other post-processed images) will show the 
enhancement pattern dynamics. . This allows radiologist 
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looking at one volume at a time to use additional 
information from all the range of dynamic volumes. The 4D 
orthogonal display provides the interface for effective 
'navigation' in 4D dynamic data. The tool can be used for 
5 both motion detection and visualisation. The ability to 
review simultaneously the subtracted image and the time- 
domain image of the whole range of subtraction images 
focused on the lesion may provide additional power to the 
interpretation process. 

10 

Modifications and variants to the above embodiments, such 
as would be readily apparent to the skilled person, are 
encompassed within the scope of the present invention. 



15 
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CLAIMS : 

1. A method of processing a time-sequence of separate 
image data sets each set comprising a plurality of image 

5 data items which each represent the location of an image 
pixel of the image subject, according to a spatially fixed 
reference frame" within which the subject is located, the 
method including the steps of: 

(a) selecting a plurality of separate image pixel 
10 locations within said reference frame; 

(b) selecting from each of a plurality of said 
separate image data sets those image data items which 
represent image pixels located at said plurality of 
separate image pixel locations; 

15 (c) generating a time-domain image according to said 

selected image data items wherein those image pixels of 
the time-domain image sharing a common image pixel 
location within said reference frame are arranged in said 
time-sequence in a line containing only those pixels. 

20 

2. A method according to Claim 1 wherein each image 
pixel is an image voxel such that each image data item- 
represents the location of an image voxel of the image 
subject according to a spatially fixed three-dimensional 

25 reference frame wherein each image data set contains 
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image data items representing image voxel locations 
spanning all three dimensions of said reference frame. 

3. A method according to any preceding claim wherein 

5 the selection of said plurality of image pixel locations 
at step (a) comprises selecting a target image pixel 
location within said reference frame and subsequently 
selecting all image pixel locations which are intersected 
by a predefined first locus intersecting said target 
10 image pixel location. 

4. A method according to Claim 3 wherein said 
predefined first locus is a line extending through said 
reference frame. 

15 

5. A method according to Claim 4 wherein said line 
extends through only one dimension of said reference 
frame and intersects all pixel locations within a column 
or a row of pixel locations within the reference frame. 

20 

6. A method according to any of preceding claims 3 to 5 
wherein said time-domain image is generated such that 
those image pixels of the time-domain image which are 
re p resen ted by image data items selected from a common 

•25 image data set are arranged in a line containing only 
those pixels. 
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7. A method according to Claim 6 in which said time- 
domain image comprises a matrix of pixels arranged to 
define columns and rows, wherein each column contains 
5 only those image pixels sharing a common pixel location 
within said reference frame, and each row contains only 
those image pixels selected from a common image data set 
or vice versa. 

10 8. A method according to any of preceding claims 3 to 7 

including: selecting a tar-get one of said plurality of 
separate image data sets and selecting therefrom those 
image data items which represent image pixels located in 
a common first plane containing said first locus; 

15 generating a first spatial-domain image according to the 
data items so selected; and, generating a composite image 
containing said first spatial-domain image and said time- 
domain image aligned such that those pixels of the first 
spatial-domain image and of the time-domain image which 

20 share a common image pixel location within said reference 
frame are arranged in-line. 

9. A method according to Claim 8 in which said 
composite image contains a second time-domain image 
25 generated according to any one of claims 3 to 7 in 

respect of a second predefined locus which intersects 



46 

said first locus at said target image pixel location 

selected at step (a) , and is substantially coplanar with 

said first locus, wherein said first spatial-domain image 

i 

and said second time-domain image are aligned such that 
5 those pixels of the first spatial-domain image and of the 
second time-domain image which share a common image pixel 
location within said reference frame are arranged in- 
line . 

10 10. A method according to Claim 9 wherein said 

first predefined locus is a line and said second 
predefined locus is a line substantially perpendicular to 
said first predefined locus. 

15 11. A method according to any of claims 8 to 10 

including: selecting from said target one of said 
plurality of separate image data sets those data items 
which represent image pixels located in a common second 
plane being non-coplanar with said first plane and 

20 containing said first locus; generating a second spatial 
domain image according to the data items so selected; 
and, generating said composite image so as to contain 
said second spatial-domain image aligned such that those 
pixels of the second spatial-domain image and of said 

25 time-domain image which share a common pixel location 
within said reference frame are arranged in-line. 
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12. A method according to any one of preceding 

claims 8 to 11 including: selecting from said target one 
of said plurality of separate image data sets those data 
items which represent image pixels located in a common 
third plane being non-coplanar with said first plane and 
containing said second locus; generating a third spatial 
domain image according to the data items so selected; 
and, generating said composite image so as to contain 
said third spatial-domain image aligned such that those 
pixels of the third spatial-domain image and of said 
second time-domain image which share a common pixel 
location within said reference frame are arranged in- 
line . 

13. A method according to Claim 12 in which said 

composite image contains a third time-domain image 
generated according to any one of claims 3 to 7 in 
respect of a third predefined locus which intersects said 
first locus and said second locus at said target image 
pixel location and is contained within one of said second 
plane and said third plane, wherein said third time- 
domain image and one of said second and said third 
spatial-domain images respectively are aligned such that 
those pixels of the third time-domain image and those 
pixels of said one of said second and said third spatial- 
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domain images which share a common image pixel location 
within said reference frame are arranged in-line. 

14. A method according to claim 12 or 13 wherein 
said first, second and third planes are all mutually 
substantially perpendicular. 

15. Apparatus for processing a time-sequence of 
separate image data sets each set comprising a plurality 
of image data items which each represent the location of 
an image pixel of the image subject according to a 
spatially fixed reference frame within which the subject 
is located, the method including the steps of: 

(a) , first selection means selecting a plurality of 
separate image pixel locations within said reference 
frame; 

(b) second selection means for selecting from each 
of a plurality of said separate image data sets those 
image data items which represent image pixels located at 
said plurality of separate image pixel locations; 

(c) image generating means for generating a time- 
domain image according to said selected image data items 
wherein those image pixels of the time-domain image 
sharing a common image pixel location within said 
reference frame are arranged in said time-sequence in a 
line containing only those pixels. 
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16. Apparatus according to Claim 15 wherein each 
image pixel is an image voxel such that each image data 
item represent the location of an image voxel of the 
image subject according to a spatially fixed three- 
dimensional reference frame wherein each image data set 
contains image data items representing image voxel 
locations spanning all three dimensions of said reference 
frame . 

17. Apparatus according to any of preceding claims 
15 and 16 wherein the first selection means is arranged 
to select a target image pixel location within said 
reference frame, and includes locus means for defining a 
first locus intersecting said target image pixel 
location, wherein the first selection means is arranged 
to select all image pixel locations which are intersected 
by said first locus. 

18. Apparatus according to Claim 17 wherein said 
predefined first locus is a line extending through said 
reference frame . 

19. Apparatus according to Claim 18 wherein said 
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* 

line extends through only one dimension of said reference 
frame and intersects all pixel locations within a column 
or a row of pixel locations within the reference frame. 

* 

5 20. Apparatus according to any of preceding claims 

17 to 19 wherein said image generating means is arranged 
to generate said time-domain image such that those image 
pixels of the time-domain image which are represented by 
image data items selected from a common image data set 

10 are arranged in a line containing only those pixels. 

21. Apparatus according to Claim 20 in which said 

image generating means is arranged to generate said time- 
domain image comprising a matrix of pixels arranged to 
15 define columns and rows, wherein each column contains 

only those image pixels sharing a common pixel location 
within said reference frame, and each row contains only 
those image pixels selected from a common image data set, 
or vice versa. 



20 
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22. Apparatus according to any of preceding claims 

17 to 21 wherein said selection means is arranged to 
select from a target one Of said plurality of separate 
image data sets those image data items which represent 
image pixels located in a common first plane containing 
said first locus; wherein the image generating means is 
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arranged to generate a first spatial-domain image 
according to the data items so selected and, to generate 
a composite image containing said first spatial-domain 
image and said time-domain image aligned such that those 
5 pixels of the first spatial-domain image and of the time- 
domain image which share a common image pixel location 
within said reference frame are arranged in-line. 

23. Apparatus according to Claim 22 wherein: said 
10 locus means is arranged to define a second locus which 

intersects said first locus at said target image pixel 
location and is substantially coplanar with said first 
locus; said second selection means is arranged to select 
all image data items the pixel locations of which are 

15 intersected by said second locus; said image' generation 
means is arranged to generate a second time-domain image 
according to image data items so selected and to generate 
a composite image which contains a second time-domain 
image in which said first spatial-domain image and said 

20 second time-domain image are aligned such that those 
pixels of the first spatial-domain image and of the 
second time-domain image which share a common image pixel 
location within said reference frame are arranged in- 
line - 

25 

24. Apparatus according to Claim 23 wherein said 
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Locus means is arranged to define said first locus as a 
line and said second locus is a line substantially 
perpendicular to said first locus. 

5 25. Apparatus according to any of claims 22 to 24 

Wherein: said second selection means is arranged to 
select from said target one of said plurality of separate 
image data sets those data items which represent image 
pixels located in a common second plane being non- 
10 coplanar with said first plane and containing said first 
locus; said image generating means is arranged to 
generate a second spatial-domain image according to the 
data items so selected, and to generate said composite 
image so as to contain said second spatial-domain image 
15 aligned such that those pixels of the second spatial- 
domain image and of said time-domain image which share a 
common pixel location within said reference frame are 
arranged in-line. 



20 26. Apparatus according to any one of preceding 

claims 22 to 25 wherein: said second selection means is 
arranged to select from said target one of said plurality 
of separate image data sets those data items which 
represent image pixels located in a common third plane 

25 being non-coplanar with said first plane and containing 
said second locus; and said image generating means is 



53 

arranged to generate a third spatial domain image 
according to the data items so selected, and to generate 
said composite image so as to contain said third spatial- 
domain image aligned such that those pixels of the third 
5 spatial-domain image and of said second time-domain image 
which share a common pixel location within said reference 
frame are arranged in-line. 

27. Apparatus according to Claim 26 wherein: said 

10 locus means is arranged to define a third locus which 

intersects said first locus and said second locus at said 
target image pixel location and is contained within one 
of said second plane and said third planes; said second 
selection means is arranged to select all image data 
15 items the pixel locations of which are intersected by 

said third locus; said image generation means is arranged 
to generate a third time-domain image according to image 
data items so selected and to generate a composite image 
which contains said third time-domain image, wherein said 
20 third time-domain image and one of said second and said 
third spatial-domain images respectively are aligned such 
that those pixels of the third time-domain image and 
those pixels of said one of said second and said third 
spatial-domain images which share a common image pixel 
location within said reference frame are arranged in- 
line. 
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28. Apparatus according to claim 26 or 27 wherein 
said second selection means is arranged to select said 

a 

first, second and third planes to be mutually 
substantially perpendicular. 

29. Apparatus according to any one of preceding 
claims 15 to 27 comprising computer means programmed to 
perform the method according to any one of claims 1 to 
13. 

30. . Computer means programmed to perform the method 
according to any one of claims 1 to 14 . 

31. A computer program product containing a 
computer program for performing the method according to 
any one of claims 1 to 14. 

32. A computer program for performing the method 
according to any one of claims 1 to 14 . 

33. An image generated according to the method of 
any one of claims 1 to 14. 



34. A method substantially as described in any 



one 
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embodiment hereinbefore with reference to the 
accompanying drawings. 

35. Apparatus substantially as -described in any 

embodiment hereinbefore with reference to the 
accompanying drawings . 




(a) 

Figure 2 
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